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ABSTRACT 



We present twenty-three transit light curves and seven occultation light curves for the ultra-short period planet WASP-43 b, in addition 
to eight new measurements of the radial velocity of the star. Thanks to this extensive data set, we improve significantly the parameters 
of the system. Notably, the largely improved precision on the stellar density (2.41 ± 0.08 p Q ) combined with constraining the age to be 
younger than a Hubble time allows us to break the degeneracy of the stellar solution mentioned in the discovery paper. The resulting 
stellar mass and size are 0.717 ± 0.025 M Q and 0.667 ± 0.01 1 R Q . Our deduced physical parameters for the planet are 2.034 ± 0.052 
M Jup and 1.036 ± 0.019 Rj„ p . Taking into account its level of irradiation, the high density of the planet favors an old age and a massive 
core. Our deduced orbital eccentricity, 0.0035 + q q^2s ' ls consistent with a fully circularized orbit. We detect the emission of the planet 
at 2.09 fim at better than 11-cr, the deduced occultation depth being 1560 ± 140 ppm. Our detection of the occultation at 1.19 /)m is 
marginal (790 ± 320 ppm) and more observations are needed to confirm it. We place a 3-cr upper limit of 850 ppm on the depth of 
the occultation at ~0.9 /im. Together, these results strongly favor a poor redistribution of the heat to the night-side of the planet, and 
marginally favor a model with no day- side temperature inversion. 

ON ■ Key words, stars: planetary systems - star: individual: WASP-43 - techniques: photometric - techniques: radial velocities 

• " 1 . Introduction Not only the origin of these planets but also their fate raises 

many questions. Their relatively large mass and small semi- 
O . There are now more than seven hundreds planets known to or- major axis imply tidal interactions with the host star that should 
■ bit around other stars than our Sun (Schneider 201 1). A signifi- lead in most cases to a slow spiral-in of the planet and to a trans- 
. . [ cant fraction of them are Jovian-type planets orbiting within 0. 1 f er of angular momentum to the star (e.g. Barker & Ogilvie 2009, 
> ■ AU of their host stars. Their very existence poses an interest- Jackson et al. 2009, Mastumura et al. 2010), the end result be- 
ing challenge for our theories of planetary formation and evo- i ng the planet's disruption at its Roche limit (Gu et al. 2003). 
lution, as such massive planets could not have formed so close The timescale of this final disruption depends mostly on the 
53 1 to their star (D'Angelo et al. 2010; Lubow & Ida 2010). Long timescale of the migration mechanism, the tidal dissipation ef- 
neglected in favor of disk-driven migration theories (Lin et al. ficiency of both bodies, and the efficiency of angular momentum 
1996), the postulate that past dynamical interactions combined losses from the system due to magnetic braking. These three pa- 
with tidal dissipation could have shaped their present orbit has rameters being poorly known, it is desirable to detect and study 
raised a lot of interest recently (e.g. Fabrycky & Tremaine 2007, in depth 'extreme' hot Jupiters, i.e. massive planets having ex- 
Naoz et al. 2011, Wu & Lithwick 2011), partially thanks to the ceptionally short semi-major axes that could be in the final stages 
discovery that a significant fraction of these planets have high of their tidal orbital decay, 
orbital obliquities that suggest past violent dynamical perturba- 
tions (e g Triaud et al 2010) The WASP transit survey detected two extreme examples of 

such ultra-short period planets, WASP- 18 b (Hellier et al. 2009) 
and WASP-19 b (Hebb et al. 2010), both having an orbital period 
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Send offprint requests to: michael.gillon@ulg.ac.be smaller than one day. Thanks to their transiting nature, these sys- 

* Based on data collected with the TRAPPIST and Euler telescopes terns were amenable for a thorough characterization that made 

at ESO La Silla Observatory, Chile, and with the VLT/HAWK-I instru- possible a study of the effects of their tidal interactions (Brown 

ment at ESO Paranal Observatory, Chile (program 086.C-0222). et al. 201 1). Notably, photometric observations of some of their 
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occultations made possible not only to probe both planets' day- 
side emission spectra but also to bring strong constraints on their 
orbital eccentricity (Anderson et al. 2010, Gibson et al. 2010, 
Nymeyer et al. 2011), a key parameter to assess their tidal his- 
tory and energy budget. The same WASP survey has recently an- 
nounced the discovery of a third ultra-short period Jovian planet 
called WASP-43 b (Hellier et al. 2011b, hereafter Hll). Its or- 
bital period is 0.81 d, the only hot Jupiter having a smaller pe- 
riod being WASP-19b (0.79 d). Furthermore, it orbits around a 
very cool K7-type dwarf that has the lowest mass among all the 
stars orbited by a hot Jupiter (0.58 + 0.05 M , T eff = 4400 ± 200 
K, Hll), except for the recently announced M0 dwarf KOI-254 
(0.59 + O.O6M , T eff = 3820 ± 90 K, Johnson et al. 2011). 
Nevertheless, Hll presented another plausible solution for the 
stellar mass that is significantly larger, 0.71 + 0.05 M . This de- 
generacy translates into a poor knowledge of the physical pa- 
rameters of the system. 

With the aim to improve the characterization of this interest- 
ing ultra-short period planet, we performed an intense ground- 
based photometric monitoring of its eclipses (transits and occul- 
tations), complemented with new measurements of the radial ve- 
locity (RV) of the star. These observations were carried out in the 
frame of a new photometric survey based on the 60cm robotic 
telescope TRAPPISlQ (77?Ansiting Planets and flanetes/mals 
.Small Telescope; Gillon et al. 2011a, Jehin et al. 2011). The 
concept of this survey is the intense high-precision photomet- 
ric monitoring of the transits of southern transiting systems, its 
goals being (i) to improve the determination of the physical and 
orbital parameters of the systems, (ii) to assess the presence of 
undetected objects in these systems through variability studies 
of the transit parameters, and (iii) to measure or put an upper 
limit on the very-near-IR thermal emission of the most highly 
irradiated planets to constrain their atmospheric properties. We 
complemented the data acquired in the frame of this program for 
the WASP-43 system by high-precision occultation time-series 
photometry gathered in the near-IR with the VLT/HAWK-I in- 
strument (Pirard et al., 2004, Casali et al. 2006) in our program 
086.C-0222. We present here the results of the analysis of this 
extensive data set. Section 2 below presents our data. Their anal- 
ysis is described in Sec. 3. We discuss the acquired results and 
drawn inferences about the WASP-43 system in Sec. 4. Finally, 
we give our conclusions in Sec. 5. 

2. Data 

2.1. TRAPPIST I+z filter transit photometry 

We observed 20 transits of WASP-43 b with TRAPPIST and its 
thermoelectric ally-cooled 2k x 2k CCD camera with a field of 
view of 22' x 22' (pixel scale = 0.65"). All the 20 transits were 
observed in an Astrodon 'I+z' filter that has a transmittance 
>90% from 750 nm to beyond 1 100 nnfl the red end of the effec- 
tive bandpass being defined by the spectral response of the CCD. 
This wide red filter minimizes the effects of limb-darkening 
and differential atmospheric extinction while maximizing stel- 
lar counts. Its effective wavelength for T e ff = 4400 + 200 K 
is A e ff = 843.5 ±1.2 nm. The mean exposure time was 20s. 
The telescope was slightly defocused to minimize pixel-to-pixel 
effects and to optimize the observational efficiency. We gener- 
ally keep the positions of the stars on the chip within a box of 
a few pixels of side to improve the photometric precision of our 



1 see http://www.ati.ulg.ac.be/TRAPPIST 

2 see http://www.astrodon.com/products/filters/near-infrared/ 



TRAPPIST time-series, thanks to a 'software guiding' system 
deriving regularly astrometric solutions on the science images 
and sending pointing corrections to the mount if needed. It could 
unfortunately not be used for WASP-43, because the star lies in 
a sky area that is not covered by the used astrometric catalogue 
(GSC1.1). This translated into slow drifts of the stars on the chip, 
the underlying cause being the imperfection of the telescope po- 
lar alignment. The amplitudes of those drifts on the total duration 
of the runs were ranging between 15 and 85 pixels in the right 
ascension direction and between 5 and 75 pixels in the declina- 
tion direction. Table 1 presents the logs of the observations. The 
first of these 20 transits was presented in HI 1. 

After a standard pre-reduction (bias, dark, flatfield correc- 
tion), the stellar fluxes were extracted from the images using the 
IRAF/DAOPHOlQ aperture photometry software (Stetson, 1987). 
For each transit, several sets of reduction parameters were tested, 
and we kept the one giving the most precise photometry for the 
stars of similar brightness as WASP-43. After a careful selection 
of reference stars, differential photometry was then obtained. 
The resulting light curves are shown in Fig. 1 and 2. 

2.2. Euler Gunn-r' filter transit photometry 

Three transits of WASP-43 b were observed in the Gunn-r' fil- 
ter (A e ff = 620.4 + 0.5 nm) with the EulerCAM CCD camera 
at the 1.2-m Euler Swiss telescope, also located at ESO La Silla 
Observatory. EulerCAM is a nitrogen-cooled 4k x 4k CCD cam- 
era with a 15' x 15' field of view (pixel scale=0.23"). Here too, 
the telescope was slightly defocused to optimize the photometric 
precision. The mean exposure time was 85s. The stars were kept 
approximately on the same pixels, thanks to a 'software guiding' 
system similar to TRAPPIST's but using the UCAC3 catalogue. 
The calibration and photometric reduction procedures were sim- 
ilar to the ones performed on the TRAPPIST data. The logs of 
these Euler observations are shown in Table 1, while the result- 
ing light curves are visible in Fig. 2. Notice that the first of these 
Euler transits was presented in HI 1. 

2.3. TRAPPIST z' filter occcultation photometry 

Five occultations of WASP-43 b were observed with TRAPPIST 
in a Sloan z' filter (A e ff - 915.9 + 0.5 nm). Their logs are pre- 
sented in Table 1. The mean exposure time was 40s. The cali- 
bration and photometric reduction of these occultation data were 
similar to the ones of the transits. Fig. 3 shows the resulting light 
curves with their best-fit models. 

2.4. VLT/HAWK-I 1.19 and 2.09 jum occultation photometry 

We observed two occultations of WASP-43 b with the cryogenic 
near-IR imager HAWK-I at the ESO Very Large Telescope in 
our program 086.C-0222. HAWK-I is composed of four Hawaii 
2RG 2048x 2048 pixels detectors (pixel scale = 0.106"), its to- 
tal field of view on the sky being 7.5' x7. 5'. We choose to ob- 
serve the occultations within the narrow band filters NB2090 
(A = 2.095 jt/m, width = 0.020 //m) and NB1190 (A = 1.186 
jum, width = 0.020 //m), respectively. The small width of these 
cosmological filters minimizes the effect of differential extinc- 
tion. Furthermore, they avoid the largest absorption and emis- 

3 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Date 


Instrument 


Filter 


N„ 


Epoch 


Baseline 
function 


<x 
[%] 


0" 120* 
[%] 






CF 


06 Dec 2010 


TRAPPIST 


I+Z 


393 


11 


pit 2 ) 


0.35, 0.36 


0.14, 0.17 


1.26, 1.30 


1.00, 1.90 


1.26, 2.47 


09 Dec 2010 


VLT/HAWK-I 


NB2090 


183 


13.5 


pit 2 ) + P H l ) + pixy') 


0.055, 0.055 


0.036, 0.036 


1.04, 1.04 


1.09, 1.14 


1.14, 1.18 


15 Dec 2010 


TRAPPIST 


I+z 


442 


22 


pit 2 ) 


0.36, 0.36 


0.12, 0.14 


0.97, 0.97 


1.03, 1.13 


1.00, 1.10 


19 Dec 2010 


TRAPPIST 


I+z 


480 


27 


pit 2 ) 


0.27, 0.28 


0.11,0.12 


1.00, 1.01 


1.00, 1.18 


1.00, 1.19 


28 Dec 2010 


TRAPPIST 


I+z 


574 


38 


P(' 2 ) 


0.32, 0.32 


0.12, 0.12 


1.01, 1.02 


1.29, 1.50 


1.30, 1.53 


28 Dec 2010 


Euler 


GunnV 


111 


38 


pit 2 ) + pixy 2 ) 


0.10, 0.10 


0.10, 0.10 


1.54, 1.64 


1.04, 1.15 


1.61, 1.90 


30 Dec 2010 


TRAPPIST 


z ' 


332 


40.5 


Pit 2 ) 


0.30, 0.30 


0.14, 0.14 


1.04, 1.04 


1.00, 1.00 


1.04, 1.04 


01 Jan 2011 


TRAPPIST 


I+z 


407 


43 


pit 2 ) + pixy 2 ) 


0.27, 0.27 


0.11,0.12 


1.00, 1.01 


1.00, 1.00 


1.00, 1.01 


06 Jan 2011 


TRAPPIST 


I+z 


273 


49 


Pit 2 ) 


0.24, 0.24 


0.11,0.11 


1.07, 1.09 


1.00, 1.38 


1.07, 1.50 


09 Jan 2011 


VLT/HAWK-I 


NB1190 


115 


51.5 


pit 2 ) + p(b 1 ) + p(xy l ) 


0.087, 0.087 


0.047, 0.048 


2.22, 2.22 


1.00, 1.00 


2.22, 2.22 


14 Jan 2011 


TRAPPIST 


I+z 


237 


59 


Pit 2 ) 


0.17, 0.18 


0.09, 0.10 


0.91,0.93 


1.04, 1.34 


0.95, 1.24 


19 Jan 2011 


TRAPPIST 


I+z 


279 


65 


Pit 2 ) 


0.18, 0.19 


0.09, 0.10 


0.91,0.95 


1.00, 2.16 


0.91, 2.05 


23 Jan 2011 


TRAPPIST 


I+z 


244 


70 


pit 2 ) + pixy 2 ) 


0.20, 0.20 


0.08, 0.09 


1.02, 1.02 


1.00, 1.00 


1.02, 1.02 


28 Jan 2011 


Euler 


Gunn-r' 


114 


76 


pit 2 ) + pixy 1 ) 


0.13, 0.13 


0.13, 0.13 


1.36, 1.42 


1.17, 1.24 


1.59, 1.76 


06 Feb 2011 


Euler 


Gunn-r' 


107 


87 


Pit 2 ) 


0.13, 0.15 


0.13, 0.15 


0.92, 1.07 


1.10, 1.86 


1.02, 1.99 


14 Feb 2011 


TRAPPIST 


I+z 


311 


97 


Pit 2 ) 


0.21,0.22 


0.11,0.12 


1.05, 1.08 


1.00, 1.42 


1.05, 1.53 


08 Mar 2011 


TRAPPIST 


I+z 


285 


124 


Pit 2 ) 


0.21,0.22 


0.11,0.13 


1.11, 1.15 


1.00, 2.00 


1.11,2.30 


10 Mar 2011 


TRAPPIST 


z' 


216 


126.5 


Pit 2 ) 


0.22, 0.22 


0.17, 0.17 


1.09, 1.09 


1.00, 1.00 


1.09, 1.09 


21 Mar 2011 


TRAPPIST 


I+z 


322 


140 


Pit 2 ) 


0.20, 0.21 


0.11,0.13 


0.83,0.88 


1.20, 2.04 


1.00, 1.80 


22 Mar 2011 


TRAPPIST 


I+z 


237 


141 


Pit 2 ) 


0.21,0.22 


0.10, 0.11 


1.10, 1.14 


1.00, 1.67 


1.11, 1.91 


28 Mar 2011 


TRAPPIST 


z' 


195 


148.5 


Pit 2 ) 


0.28, 0.28 


0.18, 0.18 


0.81,0.81 


1.03, 1.04 


0.84, 0.84 


31 Mar 2011 


TRAPPIST 


I+z 


238 


152 


Pit 2 ) 


0.24, 0.26 


0.12, 0.17 


1.04, 1.15 


1.00, 2.28 


1.04, 2.63 


02 Apr 2011 


TRAPPIST 


z' 


160 


154.5 


pit 2 )+pib<) 


0.18, 0.18 


0.11,0.11 


0.93, 0.93 


1.27, 1.30 


1.19, 1.21 


13 Apr 2011 


TRAPPIST 


I+z 


293 


168 


pit 2 ) + pixy 2 ) 


0.23, 0.23 


0.12, 0.13 


1.00, 1.01 


1.00, 1.00 


1.00, 1.01 


17 Apr 2011 


TRAPPIST 


I+z 


289 


173 


Pit 2 ) 


0.26, 0.27 


0.13, 0.15 


1.05, 1.09 


1.00, 1.26 


1.05, 1.36 


30 Apr 2011 


TRAPPIST 


I+z 


327 


189 


Pit 2 ) 


0.32, 0.33 


0.16, 0.17 


1.31, 1.33 


1.07, 1.30 


1.40, 1.73 


09 May 20 11 


TRAPPIST 


I+z 


280 


200 


Pit 2 ) 


0.19, 0.20 


0.10, 0.12 


0.98, 1.05 


1.00, 1.77 


0.98, 1.85 


1 1 May 201 1 


TRAPPIST 


z' 


214 


202.5 


Pit 2 ) 


0.21,0.21 


0.14, 0.14 


1.03, 1.03 


1.23, 1.29 


1.27, 1.33 


18 May 2011 


TRAPPIST 


I+z 


250 


211 


Pit 2 ) 


0.18, 0.18 


0.10, 0.11 


0.93, 0.96 


1.10, 1.48 


1.02, 1.42 


13Jun2011 


TRAPPIST 


I+z 


605 


243 


Pit 2 ) 


0.33, 0.34 


0.13, 0.14 


0.97, 0.99 


1.25, 1.87 


1.22, 1.84 



Table 1. WASP-43 b photometric eclipse time-series used in this work. For each light curve, this table shows the date of acquisition, 
the used instrument and filter, the number of data points, the epoch based on the transit ephemeris presented in HI 1, the selected 
baseline function (see Sec. 3.1), the standard deviation of the best-fit residuals (unbinned and binned per intervals of 2 min), and 
the deduced values for /3 W , j3 r , and CF = f3 r X f3 w (see Sec. 3.1). For the baseline function, p(e N ) denotes, respectively, a TY-order 
polynomial function of time (e = t), the logarithm of time (e = I), x and y positions (e = xy), and background (e = b). For the 
last five columns, the first and second value correspond, respectively, to the individual analysis of the light curve and to the global 
analysis of all data. 



sion bands that are present in J and K bands, reducing thus sig- 
nificantly the correlated photometric noise caused by the com- 
plex spatial and temporal variations of the background due to 
the variability of the atmosphere. As this correlated noise is the 
main precision limit for ground-based near-lR time-series pho- 
tometry, the use of these two narrow-band filters optimizes the 
photometric quality of the resulting light curves. This is espe- 
cially important in the context of the challenging measurement 
of the emission of exoplanets. 

The observation of the first occultation (NB2090 filter) took 
place on 2010 Dec 9 from 5h37 to 9h07 UT. Atmospheric con- 
ditions were very good, with a stable seeing and extinction. 
Airmass decreased from 2.1 to 1.05 during the run. Each of the 
185 exposures was composed of 17 integrations of 1.7s each (the 
minimum integration time allowed for HAWK-I). We choose to 
do not apply a jitter pattern. Indeed, the background contribution 
in the photometric aperture is small enough to ensure that the 
low-frequency variability of the background cannot cause corre- 
lated noises with an amplitude larger than a few dozens of ppm 
in our light curves, making the removal of a background im- 
age unnecessary. Furthermore, staying on the same pixels during 



the whole run allows minimizing the effects of interpixel sensi- 
tivity inhomogeneity (i.e. the imperfect flat field). The analysis 
of HAWK-I calibration frames showed us that the detector is 
nearly linear up to 10-12 kADU. The peak of the target image 
was above 10 kADU in the first images, so a slight defocus was 
applied to keep it below this level during the rest of the run. 
The mean full- width at half maximum of the stellar point-spread 
function (PSF) was 7.3 pixels = 0.77", its standard deviation for 
the whole run being 0.57 pixels = 0.06". The pointing was care- 
fully selected to avoid cosmetic defects on WASP-43 or on the 
comparison stars. 

The second occultation (NB1190 filter) was observed on 
2011 Jan 9 from 4h57 to 9h27 UT. The seeing varied strongly 
during the all run (mean value in our images =0.76", with a stan- 
dard deviation = 0.19", minimum = 0.51", maximum = 1.22") 
while the extinction was stable during the first part of the run and 
slightly variable during the second part. No defocus was applied, 
the peak of the target being in the linear part of the detector dy- 
namic in all images. Airmass decreased from 1.37 to 1.03, then 
increased to 1.13 during the run. Here too, no jitter pattern was 
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Fig. 1. Left : WASP-43 b transit photometry (12 first TRAPPIST transits) used in this work, binned per 2 min, period-folded on 
the best-fit transit ephemeris deduced from our global MCMC analysis (see Sec. 3.3), and shifted along the y-axis for clarity. The 
best-fit baseline+transit models are superimposed on the light curves. The fifth and ninth models (from the top) show some wiggles 
because of their position-dependent terms. Right: best-fit residuals for each light curve binned per interval of 2 min. 



applied. Each of the 241 exposures was composed of 17 integra- 
tions of 1.7s. 

After a standard calibration of the images (dark-subtraction, 
Hatfield correction), a cosmetic correction was applied. This cor- 
rection was done independently for each image and based on an 
automatic detection of the stars followed by a detection of out- 
lier pixels. This latter was based on a comparison of the value of 
each pixel to the median value of the eight adjacent pixels. For 
a pixel within a stellar aperture, a detection threshold of 50-<x 
was used, while it was 5-cr for the background pixels. Outlier 
pixels had their value replaced by the median value of the ad- 



jacent pixels. At this stage, aperture photometry was performed 
with DAOPHOT for the target and comparison stars, and differ- 
ential photometry was obtained for WASP-43. Table 1 presents 
the logs of our HAWK-I observations, while Fig. 3 shows the 
resulting light curves with their best-fit models. 

2.5. Euler/CORALIE spectra and radial velocities 

We gathered eight new spectroscopic measurements of WASP- 
43 with the CORALIE spectrograph mounted on Euler. The spec- 
troscopic measurements were performed between 2011 Feb 02 
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Fig. 2. Le/f : WASP-43 b transit photometry (8 last TRAPPIST transits + the 3 Eider transits) used in this work, binned per 2 min, 
period-folded on the best-fit transit ephemeris deduced from our global MCMC analysis (see Sec. 3.3), and shifted along the y-axis 
for clarity. The best-fit baseline+transit models are superimposed on the light curves (red = TRAPPIST I + z filter; green = Euler 
Gunn-r' filter). The third, ninth and tenth models (from the top) show some wiggles because of their position-dependent terms. 
Right: best-fit residuals for each light curve binned per interval of 2 min. 



and March 12, the integration time being 30 min for all of them. 
RVs were computed from the calibrated spectra by weighted 
cross-correlation (Baranne et al. 1996) with a numerical spec- 
tral template. They are shown in Table 2. We analyzed these new 
RVs globally with the RVs presented in HI 1 and with our eclipse 
photometry (see Sec. 3.3). 



3. Data analysis 



Our analysis of the WASP-43 data was divided in two steps. In 
a first step, we performed individual analyses of the 30 eclipse 
light curves, determining independently for each light curves 
the corresponding eclipse and physical parameters. The aim of 
this first step was searching for potential variability among the 
eclipse parameters (Sec. 3.2). We then performed a global anal- 
ysis of the whole data set, including the RVs, with the aim to 
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1.04 



1.04 




Fig. 3. Left : WASP-43 b occultation photometry, binned per intervals of 2 min, period-folded on the best-fit transit ephemeris 
deduced from our global MCMC analysis (see Sec. 3.3), and shifted along the y-axis for clarity. The best-fit baseline+occultation 
models are superimposed on the light curves (blue = TRAPPIST z' -filter; green = VLT/HAWK-I NB1 190 filter; red = VLT/HAWK- 
I NB2090 filter). Right: same light curves divided by their best-fit baseline models. The corresponding best-fit occultation models 
are superimposed. 



obtain the strongest constraints on the system parameters (Sec. 
3.3). 

3.1. Method and models 

Our data analysis was based on the most recent version of 
our adaptive Markov Chain Monte-Carlo (MCMC) algorithm 
(see Gillon et al. 2010 and references therein). To summa- 
rize, MCMC is a Bayesian stochastic simulation algorithm de- 
signed to deduce the posterior Probability Distribution Functions 



(PDFs) for the parameters of a given model (e.g. Gregory 2005, 
Carlin & Louis 2008). Our implementation of the algorithm 
assumes as model for the photometric time-series the eclipse 
model of Mandel & Agol (2002) multiplied by a baseline model 
aiming to represent the other astrophysical and instrumental 
mechanisms able to produce photometric variations. For the 
RVs, the model is based on Keplerian orbits added to a model 
for the stellar and instrumental variability. Our global model can 
include any number of planets, transiting or not. For the RVs 
obtained during a transit, a model of the Rossiter-McLaughlin 
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Time 

(Sy/) TZ , B -2,450,000) 


KV 

(km s" 1 ) 


(m s" 1 ) 


DC 
DO 

(kms" 1 ) 


5594.848273 


-3.921 


21 


0.021 


5604.730390 


-4.149 


19 


0.052 


5605.677373 


-3.862 


15 


0.044 


5626.663879 


-4.123 


18 


0.035 


5627.670365 


-3.759 


21 


-0.039 


5628.713064 


-3.059 


18 


0.061 


5629.686781 


-3.388 


17 


-0.009 


5632.715791 


-3.126 


22 


0.115 



Table 2. CORAL IE radial-velocity measurements for WASP-43 
(BS = bisector spans). 



effect is also available (Gimenez, 2006). Comparison between 
two models can be performed based on their Bayes factor, this 
latter being the product of their prior probability ratio multiplied 
by their marginal likelihood ratio. The marginal likelihood ra- 
tio of two given models is estimated from the difference of their 
Bayesian Information Criteria (BIC; Schwarz 1978) which are 
given by the formula: 

BIC=x 2 +klog(N) (1) 

where k is the number of free parameters of the model, N is the 
number of data points, and x 1 is the smallest chi-square found 
in the Markov chains. From the BIC derived for two models, the 
corresponding marginal likelihood ratio is given by e - ABIC / 2 , 

For each planet, the main parameters that can be randomly 
perturbed at each step of the Markov chains (called jump param- 
eters) are 

- the planet/star area ratio dF = {R p IR+) 2 , R p and /?* being 
respectively the radius of the planet and the star; 

- the occultation depth(s) (one per filter) dF occ ; 

- the parameter b' = a cos /,,//?+ which is the transit impact 
parameter in case of a circular orbit, a and i p being respec- 
tively the semi-major axis and inclination of the orbit; 

- the orbital period P; 

- the time of minimum light Tq (inferior conjunction); 

- the two parameters yje cos w an d sin w > e being the or- 
bital eccentricity and co being the argument of periastron; 

- the transit width (from first to last contact) W; 

- the parameter K 2 = K VI - e 2 P 1/3 , K being the RV orbital 
semi-amplitude; 

- the parameters Vv sin/* cos/3 and y/v sin/* sin/3, vsin/* 
and B being respectively the projected rotational velocity of 
the star and the projected angle between the stellar spin axis 
and the planet's orbital axis. 

Uniform or normal prior PDFs can be assumed for the jump and 
physical parameters of the system. Negative values are not al- 
lowed for dF, dF 0CC , V, P, T , W and K 2 . 

Two limb-darkening laws are implemented in our code, 
quadratic (two parameters) and non-linear (four parameters). 
For each photometric filter, values and error bars for the limb- 
darkening coefficients are interpolated in Claret & Bloemen's 
tables (2011) at the beginning of the analysis, basing on input 
values and error bars for the stellar effective temperature T e g, 
metallicity [Fe/H] and gravity logg. For the quadratic law, the 
two coefficients u\ and u 2 are allowed to float in the MCMC, 
using as jump parameters not these coefficients themselves but 
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Filter 




U2 


I+z 


W(0.440, 0.035 2 ) 


JV(0.180,0.025 2 ) 


Gunn-r' 


A(0.625,0.015 2 ) 


N(0.115,0.010 2 ) 



Table 3. Prior PDF used in this work for the quadratic limb- 
darkening coefficients. 



the combinations c\ — 2 x u\ + u 2 and c 2 — u\ -2xu 2 to mini- 
mize the correlation of the obtained uncertainties (Holman et al. 
2006). In this case, the theoretical values and error bars for u\ 
and u 2 deduced from Claret's tables can be used in normal prior 
PDFs. In all our analyses, we assumed a quadratic law and let u\ 
and u 2 float under the control of the normal prior PDFs deduced 
from Claret & Bloemen's tables. For the non-standard / + z fil- 
ter, the modes of the normal PDFs for u\ and u 2 were taken as 
the averages of the values interpolated from Claret's tables for 
the standard filters Ic and z' , while the errors were computed as 
the quadratic sums of the errors for these two filters. The prior 
PDFs deduced for WASP-43 are shown in Table 3. They were 
computed for T eS = 4400 + 200K, logg = 4.5 ± 0.2 and [Fe/H] 
= -0.05 ±0.17 (HI 1). 

At the first step of the MCMC, the timings of the measure- 
ments are passed to the BJD TDB time standard, following the 
recommendation of Eastman et al. (2010) that outlined that the 
commonly-used B J Dure time standard is not practical for high- 
precision timing monitoring as it drifts with the addition of one 
leap second roughly each year. 

At each step of the Markov Chains, the stellar density is 
derived from Kepler's third law and the jump parameters dF, 
b' , W, Vicosw and Vesinw (Seager & Mallen-Ornelas 2003, 
Winn 2010). Using as input values the resulting stellar density 
and values for T e ff and [Fe/H] drawn from their prior distri- 
butions, a modified version of the stellar mass calibration law 
deduced by Torres et al. (2010) from well-constrained detached 
binary systems (see Gillon et al. 201 lb for details) is used to de- 
rive the stellar mass. The stellar radius is then derived from the 
stellar density and mass. At this stage, the physical parameters 
of the planet (mass, radius, semi-major axis) are deduced from 
the jump parameters and stellar mass and radius. Alternatively, 
a value and error for the stellar mass can be imposed at the 
start of the MCMC analysis, In this case, a stellar mass value is 
drawn from the corresponding normal distribution at each step 
of the Markov Chains, allowing the code to deduce the other 
physical parameters. We preferred here to use this second op- 
tion, as WASP-43 was potentially lying at the lower edge of the 
mass range for which the calibration law of Torres et al. is valid, 
~ 0.6M o . 

If measurements for the rotational period of the star and for 
its projected rotational velocity are available, they can be used 
in addition to the stellar radius values deduced at each step of 
the MCMC to derive a posterior PDF for the inclination of the 
star (Watson et al. 2010, Gillon et al. 2011b). We did not use 
this option here despite that the rotational period of the star was 
determined from WASP photometry to be 15.6 + 0.4 days (HI 1), 
because the V sin/* measurement presented by Hll (4.0 ± 0.4 
kms -1 ) was presented by these authors as probably affected by 
a systematic error due to additional broadening of the lines. 

Several chains of 100,000 steps were performed for each 
analysis, their convergences being checked using the statisti- 
cal test of Gelman and Rubin (1992). After election of the best 
model for a given light curve, a preliminary MCMC analysis 
was performed to estimate the need to rescale the photometric 



7 



M. Gillon et al.: Thirty eclipses of the ultra-short period planet WASP-43 b 



errors. The standard deviation of the residuals was compared to 
the mean photometric errors, and the resulting ratios j3 w were 
stored. j3 w represents the under- or overestimation of the white 
noise of each measurement. On its side, the red noise present in 
the light curve (i.e. the inability of our model to represent per- 
fectly the data) was taken into account as described by Gillon et 
al. (2010), i.e. a scaling factor f3 r was determined from the stan- 
dard deviations of the binned and unbinned residuals for differ- 
ent binning intervals ranging from 5 to 120 minutes, the largest 
values being kept as /3 r . At the end, the error bars were mul- 
tiplied by the correction factor CF — /3 r x p w . For the RVs, a 
'jitter' noise could be added quadratically to the error bars after 
the election of the best model, to equal the mean error with the 
standard deviation of the best-fit model residuals. In this case, it 
was unnecessary. 

Our MCMC code can model very complex trends for the 
photometric and RV time-series, with up to 46 parameters for 
each light curve and 17 parameters for each RV time-series. Our 
strategy here was first to fit a simple orbital/eclipse model and to 
analyze the residuals to assess any correlation with the external 
parameters (PSF width, time, position on the chip, line bisector, 
etc.), then to use the Bayes factor as indicator to find the optimal 
baseline function for each time-series, i.e. the model minimiz- 
ing the number of parameters and the level of correlated noise in 
the best-fit residuals. For ground-based photometric time-series, 
we did not use a model simpler than a quadratic polynomial in 
time, as several effects (color effects, PSF variations, drift on 
the chip, etc) can distort slightly the eclipse shape and can thus 
lead to systematic errors on the deduced transit parameters. This 
is especially important to include a trend in the baseline model 
for transit light curve with no out-of-transit data before or af- 
ter the transit, and/or with a small amount of out-of-transit data. 
Having a rather small amount of out-of-transit data is quite com- 
mon for ground-based transit photometry, as the target star is 
visible under good conditions (low airmass) during a limited du- 
ration per night. In such conditions, an eclipse model can have 
enough degrees of freedom to compensate for a small-amplitude 
trend in the light curve, leading to an excellent fit but also to bi- 
ased results and overoptimistic error bars. Allowing the MCMC 
to 'twist' slightly the light curve with a quadratic polynomial in 
time compensates, at least partially, for this effect. 

For the RVs, our minimal baseline model is a scalar V r repre- 
senting the systemic velocity of the star. It is worth noticing that 
most of the baseline parameters are not jump parameters in the 
MCMC, they are deduced by least-square minimization from the 
residuals at each step of the chains, thanks to their linear nature 
in the baseline functions (Bakos et al. 2009, Gillon et al. 2010). 

3.2. Individual analysis of the eclipse time-series 

As mentioned above, we first performed an independent analysis 
of the transits and occultations aiming to elect the optimal model 
for each light curve and to assess the variability and robustness 
of the derived parameters. For all these analyses, the orbital pe- 
riod and eccentricity were kept respectively to 0.813475 days 
and zero (Hll), and the normal distribution 7V(0.58, 0.05 2 ) M 
(HI 1) was used as prior PDF for the stellar mass. For the transit 
light curves, the jump parameters were dF, b' , W and To. For the 
occultations, the only jump parameter was the occultation depth 
dF occ , the other system parameters being drawn at each step of 
the MCMC from normal distributions deduced from the values 
+ errors derived in HI 1 . 

Table 1 presents the baseline function selected for each light 
curve, the derived factors fi w , fi r , CF, and the standard deviation 



of the best-fit residuals, unbinned and binned per 120s. These 
results allow us to assess the photometric precision of the used 
instruments. The TRAPPIST data show mean values for /3„ and 
/3 r very close to 1, the I+z data having < yS w >= 1.03and<y6 r >= 
1.05, while the z' data have < fi w >= 0.98 and < /3 r >= 1.11. 
This suggests that the photometric errors of each measurement 
are well approximated by a basic error budget (photon, read- 
out, dark, background, scintillation noises), and that the level of 
correlated noise in the data is small. Furthermore, we notice that 
most TRAPPIST light curves are well modeled by the 'minimal 
model', i.e. the sum of an eclipse model and a quadratic trend in 
time. Only one TRAPPIST transit light curve requires additional 
terms in x and y, while one occultation light curve acquired when 
the moon was close to full requires a linear term in background. 
The mean photometric errors per 2 min intervals can also be 
considered as very good for a 60cm telescope monitoring a V = 
12.4 star: 0.1 1% and 0. 15% in the I+z and z' filters, respectively, 
which is similar to the mean photometric error of Euler data, 
0.12%. Euler data show also small mean values of < f3 w >= 
1.27 and < (3 r >= 1.05. Their modeling requires PSF position 
terms for 2 out of 3 eclipses, despite the good sampling of the 
PSF and the active guiding system keeping the stars nearly on 
the same pixels. This could indicate that the flatfields' quality is 
perfectible. 

For the first HAWKI light curve, taken in the NB2090 filter, 
we notice that we have to account for a 'ramp' effect (the log(f) 
term) similar to the well-documented sharp variation of the ef- 
fective gain of the Spitzer/IRAC detector at 8 //m (e.g. Knutson 
et al. 2008), and also for a dependance of the measured flux with 
the exact position of the PSF center on the chip. This position 
effect could be decreased by spreading the flux on more pixels 
(defocus), but this would also increase the background's contri- 
bution to the noise budget, potentially bringing not only more 
white noise but also some correlation of the measured counts 
with the variability of the local thermal structure. The photomet- 
ric quality reached by these NB2090 data (J3 r = 1 .09, mean error 
of 360 ppm per 2 min time interval), can be judged as excellent. 
For the NB1 190 HAWKI data, we had to include a dependance 
in the PSF width in the model. This is not surprising, consider- 
ing the large variability of the seeing during the run. We also no- 
tice for these data that our error budget underestimated strongly 
the noise of the measurements (fi w = 2.22), suggesting an un- 
accounted noise source. Still, the reached photometric quality 
remains excellent (J3 r - 1, mean error of 470 ppm per 2 min time 
interval). 

Table 4 presents for each eclipse light curve the values and 
errors deduced for the jump parameters. Several points can be 
noticed. 

- The emission of the planet is clearly (10-cr) detected at 2.09 
^m. At 1.19 /J.m, it is barely detected (~ 2.3-cr). It is not 
detected in any of the z'-band light curves. 

- The transit shape parameters b' , W, and dF, show scatters 
slightly larger than their average error, the corresponding ra- 
tio being, respectively, 1.24, 1.26, and 1.54. Furthermore, 
these parameters show significant correlation, as can be seen 
in Fig. 4. 

- Fitting a transit ephemeris by linear regression with the 
measured transit timings shown in Table 4, we ob- 
tained T(N) = 2,455,528.868289(±0.000072) + N x 
0.81347728(+0.00000060) BJD TDB . Table 4 shows (last col- 
umn) the resulting transit timing residuals (observed minus 
computed, O-C). Fig. 5 (upper panel) shows them as a func- 
tion of the transit epochs. The scatter of these timing resid- 



8 



M. Gillon et al.: Thirty eclipses of the ultra-short period planet WASP-43 b 



uals is 2.1 times larger than the mean error, 18s. No clear 
pattern is visible in the O-C diagram. 

The apparent variability of the four transit parameters is present 
in both TRAPPIST and Euler results. The correlation of the de- 
rived transit parameters is not consistent with actual variations 
of these parameters and favors biases from instrumental and/or 
astrophysical origin. This apparent variability of the transit pa- 
rameters could be explained by the variability of the star itself. 
Indeed, WASP-43 is a spotted star (Hll), and occulted spots 
(and faculae) can alter the observed transit shape and bias the 
measured transit parameters (e.g. Huber et al. 201 1, Berta et al. 
201 1). We do not detect clear spot signatures in our light curves, 
but our photometric precision could be not high enough to de- 
tect such low-amplitude structures, so we do not reject this hy- 
pothesis. On their side, unocculted spots should have a negligi- 
ble impact on the measured transit depths, considering the low 
amplitude of the rotational photometric modulation detected in 
WASP data (6+1 mmag). Still, they could alter the slope of the 
photometric baseline and make it more complex. We represent 
this baseline as an analytical function of several external param- 
eters in our MCMC simulations, but the unavoidable inaccuracy 
of the chosen baseline model can also bias the derived results. As 
described in Sec. 2.1, no pointing corrections were applied dur- 
ing the TRAPPIST runs because of a stellar catalogue problem. 
Even if the resulting drifts did not introduce clear correlations of 
the measured fluxes with PSF positions, they could have slightly 
affected the shape of the transits. These considerations reinforce 
the interest of performing global analysis of extensive data sets 
(i.e. many eclipses) in order to minimize systematic errors and 
to reach high accuracies on the derived parameters. 

3.3. Global analysis of photometry and radial velocities 

The global MCMC analysis of our data set was divided in three 
steps. For each step of the analysis, the MCMC jump parameters 
were dF, W, V , To, P, -\fecos u>, -\fesmo), three dF occ (for the 
z', NB1 190, and NB2090 filters), K 2 , and the limb-darkening co- 
efficients c\ and C2 for the I + z and Gunn-r' filters. No model for 
the Rossiter-McLaughlin was included in the global model, as no 
RV was obtained at the transit phase. For each light curve, we as- 
sumed the same baseline model than for the individual analysis. 
The assumed baseline for the RVs was a simple scalar (systemic 
velocity V y ). 

In a first step, we performed a chain of 100,000 steps with 
the aim to redetermine the scaling factors of the photometric er- 
rors. The deduced values are shown in Table 1 . It can be noticed 
that the mean < B r > for the 20 transits observed by TRAPPIST 
in / + z filter is now 1.52, for 1.05 after the individual analysis 
of the lightcurves. We notice the same tendency for the Euler 
Gunn-r' transits: < B r > goes up from 1.05 to 1.42. Considering 
the apparent variability of the transit parameters deduced from 
the individual analysis of the light curves and their correlations, 
our interpretation of this increase of the < B r > is that a part of 
the correlated noise (from astrophysical source or not) of a tran- 
sit light curve can be 'swallowed' in the transit+baseline model, 
leading to a good fit in terms of merit function but also to bi- 
ased derived parameters. By relying on the assumption that all 
transits share the same profile, the global analysis allows to bet- 
ter separate the actual transit signal from the correlated noises 
of similar frequencies, leading to larger < /3 r > values. It is also 
worth noticing that the < J3 r > for the transits, 1.50, is larger 
than the one for the occultations, 1.12, which is consistent with 
the crossing of spots by the planet during transit, but can also 
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Fig. 4. Correlation diagrams for the transit parameters deduced 
from the individual MCMC analysis of the 23 transit light 
curves. Top left: transit depth vs transit impact parameter. Top 
right: transit duration vs transit impact parameter. Bottom left: 
transit duration vs transit depth. Bottom right: TTV (observed 
minus calculated transit timing) vs transit duration. The filled 
black and open red symbols correspond, respectively, to the 
TRAPPIST and Euler light curves. 
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Fig. 5. Top: Observed minus calculated transit timings obtained 
from the individual analysis of the transit light curves as a func- 
tion of the transit epoch. The filled black and open red sym- 
bols correspond, respectively, to the TRAPPIST and Euler light 
curves. Bottom: same but deduced from the global analysis of all 
transits. 
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Epoch 


Filter 


dF 


b 




W 


To 


O-C 






[%] 




[%] 


[h] 


[BJD rDS -2450000] 


[min] 


11 


I + z 




0.682^ 


9 Q9+0.14 
7 -0.13 


1 250 +0038 

-0.034 


5537.81648!JJS 3 


-0.09 ± 0.69 


13.5 


NB2090 


0-156^ 












22 


I + z 




0-673^0 


2 50 +o.io 

J -0.11 


1 9^4+0.024 
-0.029 


5546.76490 ± 0.00020 


0.16 ±0.29 


27 


I + z 




716+ 0019 

-0.021 


2.721 ±0.077 


1.250 ±0.018 


5550.832283°°™ 


0.15 ±0.20 


38 


I + z 




645 +0 043 


9 4.-^+0 11 
^■^-0.09 


1 994+0.029 

-0.026 


5559.78038 ± 0.00021 


-0.07 ± 0.30 


38 


Gunn-r' 




0.640^ 


9 490+0.088 

—0 083 


1 101 +0.019 

—0 020 


5559.78092 ± 0.00016 


0.71 ±0.22 


40.5 


z' 


0.0903°" 












43 


I + z 




' -0.047 


2.60 ±0.17 


1 176 +0022 

1.1/ u_ 


5563.84773 ± 0.00023 


-0.12 ±0.33 


49 


I + z 




0611+ 0.042 

-0.048 


9 4 95+o.090 

-0.079 


1 210 +0020 

-0.021 


5568.728303g»» 


-0.54 ± 0.26 


51.5 


NB1190 


0.082^ 












59 


I + z 




643 +0 022 


9 504+O.O67 
^•--->-r-0.070 


1.219 ±0.029 


5576.86380 ± 0.00015 


0.51 ± 0.22 


65 


I + z 




0651 +0.025 


9 070+O.O64 
' °-0.060 


1.212 ± 0.014 


5581.744123SI 2 


-0.28 ±0.19 


70 


I + z 




0.6673™ 


2 60 +0 - 24 

-0.22 


1 999+0041 
• -0.038 


5585.8 121 l^o^oools 


0.59 ± 0.56 


76 


Gunn-r' 




0.5463°» 


9 07+0.10 
-0.09 


1 1 c 0+0.023 

1,1 °-0.020 


5590.692603K 


0.05 ± 0.25 


87 


Gunn-r' 




656 +003 ° 


9 64Q+0.089 

-0.087 


1 9^7+0.019 

-0.020 


5599.64042 ± 0.00017 


-0.57 ± 0.25 


97 


I + z 




643 +0 032 


9 45 fi +0.068 
^ JU -0 .075 


1 1 88+0.017 
1.188_ 0018 


5607.775 15!« Z\t 


-0.63 ± 0.22 


124 


I + z 




gr O+0.028 
-0.029 


2 456 +0 064 


1.213 ±0.014 


5629.740023°°°!* 


0.79 ± 0.20 


126.5 


z' 


034 +0033 

U,UJ ^_0 .023 












140 


I + z 




698 +002 ° 


9 C7f\+0.070 

072 


1.223 ±0.018 


5642.754513°°™ 


-0.86 ±0.25 


141 


I + z 




633 +0 029 


9 7 1 9+0.093 


1.240 ± 0.020 


5643.56871^S| 9 


0.18 ±0.27 


148.5 


z' 


0.0393$ 












152 


I + z 




725 +0 023 

-0.027 


2.82 ±0.11 




5652.5158l3ooo23 


-1.48 ±0.33 


154.5 


z' 


0.042^0 












168 


I + z 




606 +0 033 


2.42 ± 0.12 


1 101 +0.018 
-0016 


5665.532223°°° 22 


-0.36 ±0.32 


173 


I + z 




670 +0 029 

-0 030 


9 741 +0.098 
-0.084 


1 974+0.029 

-0.026 


5669.59946 ± 0.00025 


-0.58 ± 0.36 


189 


I + z 




596 +0 051 


9 41+0.10 
Z "^ I -0.11 


1 914+0.024 

-0.026 


5682.615923S 


0.61 ± 0.30 


200 


I + z 




677 +0 017 

"• u/ ' -0.021 


9 7QO+0.060 
JO -0.070 


1 91 0+O.OI6 


5691.56378^S! 3 


0.05 ±0.19 


202.5 


z' 


0.065^ 












211 


I + z 




0.642+™* 


2 400 +0 067 


1 9n7+°°i 4 

1 - ZU/ -0.017 


5700 51232 +0 ° 0015 

J/UU.Jl« -0.00014 


0.47 ± 0.22 


243 


I + z 




0.6873$ 


2 56 +0 19 


1 236 +005 ° 

-0.043 


5726 54436+ 00056 

J/iU.JttJU _o.00044 


1.57 ±0.81 



Table 4. Median and 1-cr errors of the posterior PDFs deduced for the jump parameters from the individual analysis of the eclipse 
light curves. For each light curve, this table shows the epoch based on the transit ephemeris presented in Hll, the filter, and the 
derived values for the occultation depth, impact parameter, transit depth, transit duration, and transit time of minimum light. The last 
column shows for the transits the difference (and its error) between the measured timing and the one deduced from the best-fitting 
transit ephemeris computed by linear regression, T(N) = 2455528.868289(0.000072) + Nx 0.81347728(0.00000060) B]D TDB . 



be explained by the much larger number of transits than occul- 
tations. 

In a second step, we used the updated error scaling factors 
and performed a MCMC of 100,000 steps to derive the marginal- 
ized PDF for the stellar density p*. This PDF was used as input 
to interpolate stellar tracks. As in Hebb et al. (2009), we used 
p* 1 ^ 3 , mostly for practical reasons. Other inputs were the stel- 
lar effective and metallicity derived by Hll, T e g = 4400 ± 200 
K and [Fe/H] = -0.05 + 0.17 dex. Those parameters were as- 
sumed distributed as a Gaussian. For each p + '^ 3 value from the 
MCMC, a random Gaussian value was drawn for stellar param- 
eters. That point was placed among the Geneva stellar evolu- 
tion tracks (Mowlavi et al. 2012) that are available in a grid fine 
enough to allow an MCMC to wander without computational 
systematic effects; a mass and age were interpolated. If the point 
fell off the tracks, it was discarded. Because of the way stars be- 



have at these small masses, within the (p+ l/3 , T e ff, [Fe/H]) space, 
the parameter space over which the tracks were interpolated cov- 
ered stellar ages well above a Hubble time. A restriction on the 
age was therefore added. Only taking ages < 12 Gyr, we obtained 
the distribution as shown in Fig. 6, refining the stellar parame- 
ters to r eff = 4520 ± 120 K and [Fe/H] = -0.01 ±0.12 dex. The 
stellar mass is estimated to M+ = 0.717 ± O.O25M . The output 
distribution is very close to Gaussian as displayed in Fig. 6. In 
addition the output p* 1/3 distribution has not been affected. Thus, 
our strongly improved precision on the stellar density combined 
with a constraint on the age (to be younger than a Hubble time) 
allows us to reject the main stellar solution proposed by Hll 
(M* = 0.58 ±O.O5M ). Because the allowed (p* 1/3 , r eff , [Fe/H]) 
space within the tracks still covers quite a large area, it is impos- 
sible at the moment to constrain the stellar age. The process is 
described in Triaud 201 1 (PhD thesis). 
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Fig. 6. Main panel: modifield Herztsprung-Russell diagram 

— 1/3 

showing the posterior PDF of p+ ' as a function of T e g after 
interpolating with the Geneva stellar evolution tracks for stars 
on the main sequence. Masses are indicated and correspond to 
the bold tracks. The solid, dashed and dotted lines correspond 

to the 1-cr, 2-cr and 3-cr contours, respectively. Secondary pan- 

—i/3 

els: marginalised PDF of p + and T e g . In both panels, there are 
three histograms: dotted: input distribution; grey: output distri- 
bution; black: output distribution for ages < 12 Gyr. 




-0.4 -0.2 0.2 0.4 
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Fig. 8. £Wer/CORALIE RV measurements period-folded on the 
best-fit transit ephemeris from the global MCMC analysis, cor- 
rected for the systemic RV, with the best-fit orbital model over- 
imposed. 



We then performed a third MCMC step, using as prior distri- 
butions for Mi,, T e ff, and [Fe/H], the normal distributions match- 
ing the results of our stellar modeling, i.e. AT(0.717, 0.025 2 ) M Q , 
/V(4520, 120 2 ) K, and M-0.01, 0.12 2 ) dex, respectively. The 
best-fit photometry and RV models are shown respectively in 
Fig. 7 and 8, while the derived parameters and 1-cr error bars are 
shown in Table 5. 



MCMC Jump parameters 
Planet/star area ratio (R,,/R+) 2 [%] 2.542+°°^ 
V = a cos i p /R* 0.656 ± 0.010 



liansit widtn W |nj 


i oaoo+0 0055 
L2089 -0.0050 


lo — 2450000 [BJD TDB \ 


cti£ C/m/; i a r\r\f\ i o 


Orbital period P [d] 


0.81347753 ± 0.00000071 


RV K 2 [ms~' d 1/3 ] 


5 1 1 

Jli - -5.0 


^[e cos oj 


020 


y/e sin uj 


-0.025 _ 0064 


dF xcrJ [ppm] 


210 -i30' < 850 (3-°") 


dF occNBli90 


790+. 32 ", < 1700 (3-cr) 


dF occNB2 090 


1560 ± 140 


ch+z 


0.983 ± 0.050 


c2 1+z 


0.065 ± 0.060 


C'l,< 


1.363 ± 0.047 




0.401 ±0.051 


Deduced stellar parameters 


ul I+z 


0.406 ± 0.026 


u2 1+z 


0.171 ±0.024 


U\f 


0.625 ± 0.024 


u2 r , 


0.112 ±0.020 


V y [km s- 1 ] 


-3.5950 ± 0.0040 


Density [p s ] 


2 4 j a+0.079 


Surface gravity log g+ [cgs] 


4645+ o.on 

t.utj 0()10 


Mass M* [M Q ] 


0.717 ±0.025 


Radius [R e ] 


667 +0 01 ° 

"•""'-0.011 


Teff [K]» 


4520 ± 120 


[Fe/H] [dex]" 


-0.01 ±012 


Deduced planet parameters 


RV K[ms-' ] 


547.9+^ 


R P /R* 


1 ->QzL<; +000076 
"• 1J7 ^ -0.00077 


K [R t ] 


0.658O$J» 


b oc [R*] 


655 +0 012 


T oc - 2450000 [BJD TDB ] 


5726 95069+ 000084 

—0.00078 


Orbital semi-major axis a [AU] 


0.01526 ± 0.00018 


Roche limit a R [AU] 


0.00768 ± 0.00016 


a/a R 


1 986 +003 ° 


a/R* 


4 918 +0053 

^• 7io -0.051 


Orbital inclination L [deg] 


82.33 ± 0.20 


Orbital eccentricity e 


0.0035^;™g, < 0.0298 (3-cr) 


Argument of periastron a> [deg] 


_ 32 +ii5 


Equilibrium temperature T eq [K]* 


1440!* 


Density p p [p Jup ] 


1 S9fi+ 0084 
1.8ZO _ 007g 


Density p p [g/cm 3 ] 


1 ^77+0.063 
1 " J ' '-0.059 


Surface gravity log g p [cgs] 


3 672 +0013 
j.u/z, 0012 


Mass M p [M Jlip ] 


2 034 +0 052 

005 j 


Radius R p [R Sup ] 


1.036 ±0.019 



Table 5. Median and 1-cr limits of the posterior marginalized 
PDFs obtained for the WASP-43 system derived from our global 
MCMC analysis. "From stellar evolution modeling (Sect. 3.3). 
fc Assuming A-0 and F=\ . 
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Fig. 7. Left: TRAPPIST I + z (top) and Euler Gunn-r' (bottom) transit photometry period-folded on the best-fit transit ephemeris 
from the global MCMC analysis, corrected for the baseline and binned per 2 min intervals, with the best-fit transit models over- 
imposed. Right: same for the occultation photometry obtained by TRAPPIST in z' filter (top), and by VLT/HAWK-I in NB1190 
(middle) and NB2090 (bottom) filters, except that the data points are binned per 5 min intervals. 



Global analysis of the 23 transits with free timings 

As a complement to our global analysis of the whole data set, we 
performed a global analysis of the 23 transit light curves alone. 
The goal here was to benefit from the strong constraint brought 
on the transit shape by the 23 transits to derive more accurate 
transit timings and to assess the periodicity of the transit. In this 
analysis, we kept fixed the parameters Tq and P to the values 
shown in Table 5, and we added a timing offset as jump param- 
eter for each transit. The orbit was assumed to be circular. The 
other jump parameters were dF, W, b', and the limb-darkening 
coefficients c\ and C2 for both filters. The resulting transit tim- 
ings and their errors are shown in Table 6. Fitting a transit 
ephemeris by linear regression with these new derived transit 
timings, we obtained T(N) = 2,455,528.868227(+0.000078) + 
jVx0.81347764(±0.00000065)BJD rz , B . Table 6 shows (last col- 
umn) the resulting transit timing residuals (observed minus com- 
puted, O-C). Fig. 5 (lower panel) shows them as a function of the 
transit epochs. The scatter of these timing residuals is now 1.8 
times larger than the mean error, 19s. Comparing the errors on 
the timings derived from individual and global analysis of the 
transit photometry, one can notice that the better constraint on 
the transit shape for the global analysis improves the error of a 
few timings, but that most of them have a slightly larger error 
because of the better separation of the actual transit signal from 
the red noise (see above). 



4. Discussion 

4. 1 . The physical parameters of WASP-43 b 

Comparing our results for the WASP-43 system with the ones 
presented in Hll, we notice that our derived parameters agree 
well with the second solution mentioned in H 1 1 , while being sig- 
nificantly more precise. WASP-43 b is thus a Jupiter-size planet, 



twice more massive than our Jupiter, orbiting at only ~0.015 
AU from a 0.72 M main-sequence K-dwarf. Despite having the 
smallest orbital distance among the hot Jupiters, WASP-43 b is 
far from being the most irradiated known exoplanet, because of 
the small size and temperature of its host star. Assuming a solar- 
twin host star, its incident flux ~9.6 10 8 erg s cirT 2 would cor- 
respond to a P = 2.63 d orbit (a = 0.0374 AU), i.e. to a rather 
typical hot Jupiter. With a radius of 1.04 R Jup , WASP-43 b lies 
toward the bottom of the envelope described by the published 
planets in the R p vs. incident flux plane (Fig. 9). Taking into 
account its level of irradiation, the high density of the planet 
favors an old age and a massive core under the planetary struc- 
ture models of Fortney et al. (2010). Our results are consistent 
with a circular orbit, and we can put a 3-cr upper limit <0.03 
to the orbital eccentricity. Despite its very short period, WASP- 
43 b has a semi-major axis approximatively twice larger than its 
Roche limit a«, which is typical for hot Jupiters (Ford & Rasio 
2006, Mastumura et al. 2010). As noted by these authors, the 
inner edge of the distribution of most hot Jupiters at ~2 cir fa- 
vors migrational mechanisms based on the scattering of planets 
on much wider orbit and the subsequent tidal shortening and cir- 
cularization of their orbits. Thus, WASP-43 b does not appear to 
be a planet exceptionally close to its final tidal disruption, unlike 
WASP-19b that orbits at only 1.2 a R (Hellier et al. 201 la). 



4.2. The atmospheric properties of WASP-43 b 

We have modeled the atmosphere of WASP-43 b using the meth- 
ods described in Fortney et al. (2005, 2008) and Fortney & 
Marley (2007). In Fig. 10, we compare the planet-to-star flux 
ratio data to three atmosphere models. The coldest model (or- 
ange) uses a dayside incident flux decreased by 1/2 to simulate 
the loss of half of the absorbed flux to the night side of the planet. 
Clearly the planet is much warmer than this model. In blue and 
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0.24 

:0.32 
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0.32 
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0.27 

:0.22 

0.40 
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= 0.33 

0.33 

0.33 
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Table 6. Median and l-cr errors of the posterior PDFs de- 
duced for the timings of the transits from their global analy- 
sis. The last column shows the difference (and its error) be- 
tween the measured timing and the one deduced from the best- 
fitting transit ephemeris computed by linear regression, T(N) = 
2455528.868227(0.000078) + N x 0.81347764(0.00000065) 



red are two models where the dayside incident flux is increased 
by a factor of 4/3 to simulate zero redistribution of absorbed flux 
(see, e.g., Hansen 2008). The red model features a dayside tem- 
perature inversion due to the strong optical opacity of TiO and 
VO gases (Hubeny et al. 2003, Fortney et al. 2006). The blue 
model is run in the same manner as the red model, but TiO and 
VO opacity are removed. 

The blue and red models are constructed to maximize the 
emission from the dayside of the planet. Figure 10 shows that 
only such bright models could credibly match the data points. 
While a dayside with no temperature inversion is slightly favored 
by the 1 .2 pm data point, it is difficult to come to a firm conclu- 
sion. Day-side emission measurements from Warm Spitzer will 
help to better constrain the atmosphere as well. They will also 
allow getting an even smaller upper-limit on the orbital eccen- 
tricity. 

Fortney et al. model fits to near infrared photometry of other 
transiting planets (e.g., Croll et al. 2010) have generally fa- 
vored inefficient temperature homogenization between the day 
and night hemispheres, although the warm dayside of WASP- 
43 b appears to be at the most inefficient end of this contin- 
uum. The apparent effeciency of temperature homogenization is 
expected to vary with wavelength. In particular, near infrared 




10 9 10 s 10 7 

Incident flux [erg s 1 cm" 2 ] 



10° 



Fig. 9. Planetary radii as a function of incident flux. WASP- 
43 b is shown as a black square. Grey filled circles are 
Kepler planetary candidates (see Demory & Seager 2011). 
Transiting giant planets previously published, and mostly from 
ground-based surveys, are shown as grey triangles. The rel- 
evant parameters R p , R s , T e jf and a have been drawn from 
http://www.inscience.ch/transits on August 29, 201 1. 
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Fig. 10. Model planet-to-star flux ratios compare to the three 
data points. The data are green diamonds with l-cr error bars 
shown. The orange model assumes planet-wide redistribution of 
absorbed flux. The red and blue models assume no redistribution 
of absorbed flux, to maximize the day-side temperature. The red 
model includes gaseous TiO and VO and has a temperature in- 
version. [See the figure inset.] The blue and orange models have 
TiO and VO opacity removed, and do not have a temperature 
inversion. For each model, filled circles are model fluxes aver- 
aged over each bandpass. The data slightly favor a model with 
no day-side temperature inversion. 



bands, at minima in water vapor opacity, are generally expected 
to probe deeper into atmospheres than the Spitzer bandpasses. 



4.3. Transit timings 

Dynamical constraints can be placed on short orbits companion 
planets from the 23 transits obtained in this study. The linear fit 
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Fig. 11. Detectivity domain for a putative WASP-43 c planet, as- 
suming e c = (black) and e c = 0.05 (red). The solid curves de- 
limit the mass-period region where planets yield maximum TTV 
on WASP-43 b above 114 s (3-cr detection based on the present 
data). The dotted curves show the 1-cr threshold. Nearly hori- 
zontal solid, dashed and dotted lines shows RV detection limits 
forRV semi-amplitude K-5, 10 and 15 ms~' respectively. 

to the transit timings described above yields a reduced ;f 2 of 4.6 
and the rms of its residuals is 38s. Two transits (epochs 140 and 
152) have a O-C different from zero at the ~4-cr level. The most 
plausible explanation for the significant scatter observed in the 
transit timings is systematic errors on the derived timings due to 
the influence of correlated noise. Another potential explanation 
is asymmetries in the transit light curves caused by the crossing 
of one or several star spots by the planet. In such cases, the fitted 
transit profile is shifted with time, producing timing variations. 

We also explored the detectability domain of a second planet 
in the WASP-43 system. To this end, we followed Agol et 
al. (2005) and used the Mercury n-body integrator package 
(Chambers 1999). We simulated 3-body systems including a sec- 
ond companion with orbital periods ranging between 1.3 and 
50 days, masses from 0.1 M e to 2.0 Mj up and an orbital ec- 
centricities of e c — and e c = 0.05. For each simulation, we 
computed the rms of the computed transits timings variations. 
The results are shown on Fig. 11 where the computed 1-cr (red) 
and 3-cr (black) detection thresholds are plotted for each point in 
the mass-period plane. We also added RV detectability threshold 
curves based on 5 (solid), 10 (dash) and 15 (dot) ms~' semi- 
amplitudes (K). 

Based on the present set of timings, the timings variations 
caused by a 5 Earth mass companion in 2:1 resonance would 
have been detected with 3-cr confidence, while unseen with ra- 
dial velocities alone. 

5. Conclusions 

In this work we have presented 23 transit light curves and 7 oc- 
cultation light curves for the ultra-short period planet WASP- 
43 b. We have also presented 8 new measurements of the ra- 
dial velocity of the star. Thanks to this extensive data set, 
we have significantly improved the parameters of the system. 
Notably, our strongly improved precision on the stellar den- 
sity (2.41 + O.O8p ) combined with a very reasonable constraint 
on its age (to be younger than a Hubble time) allowed us to 
break the degeneracy of the stellar solution mentioned by HI 1. 
The resulting stellar mass and size are 0.717 + 0.025M o and 



0.667 + 0.01 IR Q . Our deduced physical parameters for the planet 
are 2.034 ± 0.052M Jup and 1.036 ± 0M9R Jup . Taking into ac- 
count its level of irradiation, the high density of the planet favors 
an old age and a massive core. Our deduced orbital eccentricity, 
0.0035*qq959, is consistent with a fully circularized orbit. 

The parameters deduced from the individual analysis of the 
23 transit light curves show some extra scatter that we attribute 
to the correlated noise of our data and, possibly, to the crossing 
of spots during some transits. This conclusion is based on the 
correlation observed among the transit parameters. These results 
reinforce the interest of performing global analysis of extensive 
data sets in order to minimize systematic errors and to reach high 
accuracies on the derived parameters. 

We detected the emission of the planet at 2.09 /mi at better 
than 1 1-cr, the deduced occultation depth being 1560+ 140 ppm. 
Our detection of the occultation at 1.19 fim is marginal (790 ± 
320 ppm) and more observations would be needed to confirm 
it. We place a 3-cr upper limit of 850 ppm on the depth of the 
occultation at ~0.9 /urn. Together, these results strongly favor a 
poor redistribution of the heat from the dayside to the nightside 
of the planet, and marginally favor a model with no day-side 
temperature inversion. 
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